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PREFACE 


This report describes changes made to the NASTRAN Thermal Analyzer 
for Variance Analysis and plotting of HBDY Elements. The theory and 
the User Information are given. Three demonstration problems were used 
to test the capability. 



IOTRODilJCTION 


The NASTRAN Thermal Analyzer has been extended to do variance analysis 
and plot the thermal boundary elements. The capability to compute the 
temperatures resulting from specified flux, boundary temperatures, radia- 
tion, etc., had been installed previously^. 

The objective of the variance analysis addition is to assess the 
sensitivity of temperature variances resulting from uncertainties inherent 
in input parameters for heat conduction analysis. The plotting capability 
provides the ability to check the geometry (location, size and orientation) 
of the boundary elements of a model in relation to the conduction elements. 

VARIANCE ANALYSIS 
Method 

Variance analysis is the study of uncertainties of the computed 
results as a function of uncertainties of the input data. To study this 
problem using NASTRAN, a solution is made for both the expected values of 
all inputs, plus another solution for each uncertain variable. A variance 
analysis module subtracts the results to form derivatives, and then can 
determine the expected deviations of output quantities.. 

For heat transfer analysis, any real number in the Bulk Data Deck 
may be treated as an uncertain item. This can include conductivity, heat 
capacity, film coefficients, emissivities , applied heat, boundary tempera- 
tures, etc. The output quantities can include temperatures, heat of con- 
straint at specified temperature boundaries , temperature gradients and 
fluxes v/ithin heat conduction elements, and radiation exchange between 
boundary elements. 
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Theory and Definitions, Variance Analysis 


Variance analysis is a study of the statistical properties of random 

2 

variables. Some terms will be defined using a notation similar to Feller . 
Let U be a random variable, and E(U) be the expectation (or expected 
value) of the random variable. The mean value of a random variable U is 
given by 

\x = E(U) (1) 

A new random variable can be defined, U-p, which is the deviation from 
the mean . The expectation of the square of the deviation from the mean 
is called the variance. 


Var(U) = E((U-p) 2 ) (2) 

The positive square root of the variance is called the standard deviation . 

a = (Var (U) )^ (3) 

The standard deviation is a measure of the scatter of the values taken by 
the random variable, and can also be interpreted as the width of the 
probability distribution curve. Other measures of scatter are possible 
(which differ from the standard deviation by a constant factor). The 
term deviation will be used for any measure of scatter; for example, the 
range which contains ninety-nine percent of the data points. The deviation 
will be written Dev(U). 

Variance analysis applied to structural analysis deals with the 
deviations of computed quantities as a result of deviations of data items. 
In the general case, there may be several data items which have uncertainty 
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and for which a mean value and deviation are known. The question is to . 
determine the mean and deviation of the output quantities. Let the data 
items be D^, D 2 / ..., D^. Let the output quantity be U. From a standard 
structural analysis, the result is given by 

U = U(D X , D 2 , .... D n ) (4) 


If the data items are independent, or uncovvetated , then the deviation of 
the result is given by 


Dev(U) = 



Dev (D . )\ 


(5) 


The basic requirement for Equation 5 to hold is that the deviations of 
the. input quantities must be sufficiently small that the output quantities 
are linear functions of the inputs. 


The derivatives will be computed from 


3U 

3D. 



( 6 ) 


-where AU is the computed result with the jth quantity changed minus the 
basic result, 6 is a user specified parameter. The user must select the 
size of the change of the input parameter sufficiently large that Equation 
6 is evaluated accurately. Usually a one to ten percent change is about 
the right size. If the change is too small,, accuracy will be lost in the 
subtraction needed to compute the numerator. If the change is too large, 
nonlinear behavior will be seen. 
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Plotting HBDY Elements 

The boundary condition elements for the NASTRAN heat conduction 

analyzer can be plotted for undeformed structure plots. This has been 

integrated into the structural plotter. General user instructions are 

3 . 

given in the NASTRAN User’s Manual , Section 4. Some special instruc- 
tions are given in the MSC/NASTRAN Application Manual^, Sections 2.3 and 
2.5. This section only describes changes from the documentation refer- 
enced above. 

There are several types of HBDY elements, as follows: 


Type 


Number of 
Primary 
Grid Points 


Normals 

Available 


P0INT 

1 

yes 

LINE 

2 

yes 

AREA3 

3 

yes 

AREA4 

4 

yes 

REV 

2 

no 

ELCYL 

2 

no 


The secondary grid points are used for ’’ambient” conditions, and are 
ignored by the plotter. Type P0INT must have a nonzero associated area 
(see AF on the associated PHBDY data card) and a defined normal direction 
(see VI, V2, V3 on the CHBDY data card) to be plotted. It is plotted as 
a hexagon with approximately the correct area. Type LINE must have a non- 
zero width (see AF on the associated PHBDY data card) and a normal defined 
in order to be able to plot. The others plot via their associated geo- 
metry defined by primary grid points. 
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BULK DATA DECK • 


Input Date'. Cura 

Description: Defines functional relationships for bulk data changes. Used for 

i variance analysis. 




1 VARY iAjrffM 

a 4 

Fi'eLd [K 

A 

■■ I 1 

e 

ID1 

ID2 

ID3 


1 VARY TmAT4 

3 jl_ 



1.0 

7 

- ■ ■ ■ 

THRU 

_ _ 

9 

— — 


fjeld 

ittol 

FIELD 
K ' 

A, B 

IDl thru ID? 


Contents 

f/ame of bulk data cards. to be modified (BCD). 

Location of data items to be changed (integer > 0). 

Parameter index (integer >0). < 

Coefficients (real). « - 

Identifiers of data cards to be modified (integer) or ID2 may 
be HTHRlir See Remark 2. 


Remarks: 1. If ID >0, the i dentifiers ID1 thru ID3 refer to field 2 of the card type 

'■ referenced by MU EM. If field 2 is not unique, then the following method 

is used. If ID<0, then ID is the sorted parent card number within the 
MNEM to change. The first card of a MNEM is numbered 1, and continuation 
cards are not 03 tX/jted. For a given hV'/EM, positive ano negative io's can 
not be mixed. 

2. Several 1VARY cards may reference the same MNEM and FIELD, but they must 
have unique values of K. If, for example, three cards appeared with fields 
2, 3, 7, 8, and 9 as above and fields 4, 5, and 6 values of 


K 

A 

B 

1 

.01 

1.0 

2 

.02 

- 0. 

3 

.04 

1.0 


then the function defined would be D = [D Q ( 1 + . 01 S , ) + ,02S^] (1+. 04S^ ) . 

Here, D refers to the item in field 3 of the MAT 4 data cards whose ID's 
are 7, 8, and 9. D Q is the value found on the card. The S's are variance 

parameters . In the general case, a sequence is- defined by 


'°k-A 1+a kV 

• 

c r 

O 

D k-1 + a k S k 

b k = o , 


u k 

and then D = last value of in sequence. 

Figure 1. IV ARY Bulk Data Card 


- 5 - 



For each case, it is required that the data item be expressed in terms of 
a parameter, S whose mean is zero and deviation, is unknown. 


k = . 167 (1+sp , 

Dev(Sp = .025 

(9) 

= 350. +S 2 , 

Dev(S 2 ) = 20. 

(10) 


The above two examples are expressed in terms of percentage and absolute 
deviations. The 1VARY Bulk Data card (see Figure 1) defines the functional 
relationship between data items and parameters. 

To illustrate the format of the data card, let us assume that the 
basic value of k is in field 3 of a MAT 4 data card of ID - 201; and that 
the basic value of T^ is given in fields 4, 6 and 8 of a TEMP Bulk Data 
card that is the first TS4P card in the deck. The functions are given by 


1VARY 

KNEW. 

FIELD K 

A B 

ID1 

ID2 IDS 

IV ARY 

MAT 4 

3 1 

1.0 1.0 

201 


1VARY 

TEMP 

4 2 

1.0 

-1 


1VARY 

TEMP 

6 2 

1.0 

-1 


1VARY 

TEMP 

8 2 

1.0 

-1 


Examination of 

the data 

card description shows 

that any 

data item can 

depend on more 

than one 

parameter. 

Another important feature is that 

many items may depend upon a single parameter. 

such as 

for the tempera- 

tures at three points in the above example. If (for example) the three 


temperatures were independent, then separate values of k should be 
assigned for the temperatures, since they depend upon independent para- 
meters. 
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BULK DATA DECK 


Input Data Card 1 VAC/ 

Description: Defines functional relationships for bulk data changes. Used for 

variance analysis. 


For mat an d Example: 


1 

flVARY 

C O 

l-'ilEM | FIELD 

K 

A 

B 

IDl 

ID2 

ID3 

1 

* 1VARY 

MAT 4 |3 

1 

.01 

1.0 

7 

THRU 

9 

1 


I! Field 
mt M 
FIELD 
K 

A, B 

ID1 thru ID3 


Contents 

Name of bulk data cards. to be modified (BCD). 

Location of data items to be changed (integer >0). 

Parameter index (integer >0). 

Coefficients (real). 

Identifiers of data cards to be modified (integer) or ID2 may 
be “THRU.** See Remark 2. 


Remarks: 1. If ID > 0, the identifiers ID1 thru ID3 refer to field 2 of the card type 

referenced by MNEM. If field 2 is not unique, then the following method 

is used. If ID < 0, then ID is the sorted parent card number within the 
MNEM to change. The first card of a MNEM is numbered 1, and continuation 
cards are not counted. For a given MNEM, positive and negative ID ! s can 
not be mixed. 

2. Several 1VARY cards may reference the same MNEM and FIELD, but they must 
have unique values of K. If, for example, three cards appeared with fields 
2, 3, 7, 8, and 9 as above and fields 4, 5, and 6 values of 


K 

A 

B 

1 

.01 

1.0 

2 

.02 

. 0 . 

3 

.04 

1.0 


then the function defined would be D = [D o (l + .01S-j ) + .C^SglO + .CMS^) . 

Here, D refers to the item in field 3 of. the MAT 4 data cards whose -ID's 
are 7, 8, and 9. D Q is the value found on the card. The S's are variance 

parameters. In the general case, a sequence is defined by 


D k-i< 1+a kV v° 


D k-l + a k s k b k = 0 ’ 


I and then D = last value of D k in sequence. 

i 

Figure 1. IV ARY Bulk Data Card 


i 
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the data cards would be 


IP ARM JRUN K1 SKI K2 SK2 K3 SK3 

IP ARM 1 1 .01 

1PARM 2 2 .01 

This data card allows several parameters to be varied for any JRUN. 
Thus, one could simultaneously increase the conductivity by twenty 
percent and the boundary temperature by one hundred degrees by using 
the 1PARM data card. For variance analysis, the parameters are usually 
varied one at a time. 


The deviations of the parameters are supplied on a VARIAN Bulk Data 
card (Figure 3), and a scale factor, 6, is supplied on a PARAM Bulk Data 
card with the name DELTA. Let U be any output quantity and be the 
data items. These are related by Equation 5. NASTRAN computes 


Dev (U) 



(13) 


The value for Dev(S^) is determined as follows: 


Dev(S k ) = g- 


Dev (D k ) 



(14) 


Frequently, one can choose all S k (on IP ARM data cards) to be equal to 6, 
making the first factor 1. The values for the two special cases are 
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BULK DATA DECK 


Input Data Card 1 PARM 

Description : Defines parameter values for modifying bulk data. Is used with 1 VARY 

~~~ data cards. 

format and Example: 


2 . 3 - 4 5 6 7 9 , JO 


l 

1 FARM 

c. 

JRUN 

Kl 

i 

SKI 

1/ 9 

fi\ C 

SK2 

K3 

SK3 

>< 


IP ARM 

3 

1 

.97 

3 

1.0 



X 

, — 


Field • Contents 


JRUN The run index for which the parameters are to be used 

(integer > 0) . 

Kl, K2, K3 • Index for nonzero parameters (integer). 

SKI , SK2, SK3 ' The values of the parameters (real). 

. Remarks : 1. The default value is given by 

$ fl.O k = ORUN 

k (0.0 k f JRUN . 

These values are used if no 1PARM cards are found for any value JRUN 

2 . If there are more than three nonzero S^‘s for any va^ue of JRUN, additional 
1PARM data cards may be introduced with the same value for JRUN. 


Figure 2. IP ARM Bulk Data Card 



BULK DATA CARD 


Input Data Card: VARTAN • • 

Description: Defines the values of the deviation of the parameters, for use in 

variance analysis. 



F ield 

DEL SI, DEL S2 , etc. 


Contents 

The deviation of the parameters used in variance analysis. 
(Real, a blank field ends the card) 


Remarks: 1. The deviation of any output quantity 0 will be computed by 


DEVIATION (0) = E(||- AS k ) 

k 


where 30/SS k are partial derivatives and AS k are values supplied on tins 
data card. 

2. Deviations will be computed only if JRUN (the loop index) is equal to the 
number of entries on this card. Only one logical card may be in the deck. 


Figure 3. VARIAN Bulk Data Card 
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The value for Dev(S^) is determined as follows: 


Dev(S k ) = —• 


Dev (D k ) 



( 14 ) 


Frequently, one can choose all $ k (on 1VARY data cards) to be equal to 6, 
making the first factor 1. The values for the two special cases are 


Dev (S k ) 



Dev(D k ) 

W 


Dev(D k ) 


percentage 

ivrm 


absolute 

form 


( 15 ) 


For the above example of data (see Equations 9 and 10), the Bulk Data 
cards used are 


PARAM DELTA .01 

VARIAN .025 20. 

All of the above Bulk Data cards must be in the Bulk Data Deck on 
restart, either as data supplied to the checkpoint run or as input data 
to the restart. The user is responsible to assure that data variations 
are large enough that numerical differences are not negligible, and yet 
small enough that the output varies linearly. 

A special parameter MESH can be given the value YES to cause printout 
of VIEW data. This uses a special VIEW Bulk Data card defined by GSFC. 
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The Case Control Deck 


Standard output requests are made. The values of the output for the 
changed parameters will be output. Wherever possible (i.e., when output 
was requested for the original data and the changed data) the derivatives 
will be output. The word "DERIVATIVES” and the value of the index JRUN 
are automatically printed to identify derivatives. On the final execution 
derivations are computed which are labeled VARIANCES." 

Standard plot requests can be made which will plot the HBDY elements. 
A new PL0T parameter "N0RMALS" is used to request vectors normal to the 
HBDY elements. The mneumonic HBDY will be recognized in set requests, 
and HB is used to label elements in plots. 

The Executive Control Deck 

A special ALTER is required to include the VARIAN and PLTHBDY modules 
These ALTERS are supplied for the three heat rigid formats. The user 
should request DIAG 40 in order to get a message about data cards affected 
by 1VARY and IP ARM data cards. For JRUN = 0, request CHKPNT YES. 

The NASTRAN Card • 

The NASTRAN card can be used to control looping. The initial and 
final values of JRUN are specified as SYSTEM(61) and SYSTEM(62). Using 
looping (called the automatic mode), the only data items that can be 
changed are those specified by 1VARY and 1PARM. Single executions can be 
made with JRUN specified by SYSTEM (61) where any data item can be changed. 
For the initial checkpoint run, JRUN = 0, so no SYSTEM(61) need be 
supplied (i.e., the default is 0). 
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The Job Submission 


It is suggested that the MERGE3 program be used with NASTRAN. This 

j 

will allow the following: 

1. $VIEW data cards may be in the input data deck. MERGE3 will 
change these to VIEW since NASTRAN ignores all DATA cards with 
$ in column 1. 

2. Rigid Format ALTERS can be merged into the deck. 

3. Checkpoint Dictionaries can be merged into the deck. 

Details of how to use MERGE 3 are given in the MSC/NASTRAN Application 
Manual . 

During the first run (JRUN = 0), the user must request tapes NPTP 
and INPT. Punch output of a checkpoint dictionary will occur. For 
later runs, the restart tape (now called 0PTP for read only), the tape 
INPT (for read and write), and the checkpoint dictionary must be supplied. 

DEMONSTRATION PROBLEMS 

Three demonstration problems are used to illustrate variance analysis 
and plotting of HBDY elements. The problems were chosen to incorporate 
many features of heat conduction problems. The actual runs have not been 
incorporated in this report. 

Variance Analysis, Nonlinear, Steady 

The model is shown in Figure 4. The variance values are given in 
Figure 5. Two computer runs are required. 

Run 1. This is a CHKPNT run using the basic values. An ALTER for 
Heat Sol 3 is used, which saves solutions on a user tape INPT for 
variance calculations. 
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4 - 


A bent plate of a uniform thickness (6) is composed of two 
materials A and B. The system is subjected to thermal 
boundary conditions consisting - of a directional energy flux 
{Q V ect) ' an internal energy source (Q V ol ) • a prescribed 
temperature boundary (T b ) , and a convective coupling with a 
fluid of constant temperature (Tf) . Diffuse-gray radiations 
.among surface elements and to an environment of temperature T<^> 
take place on the upper side of the plate while the underside 
is insulated perfectly. All dimensions, thermophys ical 
properties, values of various parameters including . the radiation 
exchange factors (A±Fij) of the finite-element model are given. 
Determine partial derivatives of the grid point temperatures 
with respect to each uncertain parameter and the resulting 
temperature intervals of each 'temperature based on 


(1) the specified variances of parameters with entries 
directly in absolute values or in percentage as given in 
following data list; 


(2) the uncertainty effect on temperatures of each parameter (s) 
to be evaluated for perturbations (As/S) of 5 %. for every 
parameter excluding Fj_j. _ 


Data: 

A • B 


k Watt/cm-k. ,167 (+2.5%) 

pep watt-sec/cm^-k 3.966 ( + 0.15) 

e IR = a IR * 9 ■ ( ±3% ) 

a solar - 95 ' 


3.933 • (±0.04) 
3.449 
. 8 ( ±3%) 

.75 


h s 2.5 x 10“ 2 watt/cm 2 -k (+15%) 
Q V ect = 0*13 watt/cm 2 (+5%) 

Q V ol = 0.52 watt/cm^ 


T b a 350 °K (±20) 

T f = 1200 °K (±30) 
Tabs' = 273.00 °K 


Af Ffj given in the listing of RADMTX (±7% to those 
elements only that i — 5 , 6, 15, 16 and 
j = 7, 8 , 17, 18) 

C T= 5.6696 x 10“ 12 watt/cm 2 -K 4 . 


Figure 5. V^-t^ice Values fur Demonstration Problems 1 and 3 
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Run 2 . This is a restart, demonstrating two loops through variance 
analysis. For additional loops, it would be necessary to change the 
"NASTRAN" card and the VAR IAN Bulk Data card. The ALTER for Heat Sol 
3 is used, along with special Bulk Data cards for variance analysis. 
For this example, the special cards are the input data deck; however, 
they could have been in the data deck of the original run. A user 
message is supplied (request DIAG 40 in the Executive Control Deck) 
which tells which items have been changed. This message appears 
after the sorted data echo. Restart bits are set which then deter- 
mine the modules which must be executed to form the changed solution. 
Output requests produce values for the new solution, the derivatives 
of the solution, and the deviations. 

Plotting HBDY Elements 

A demonstration problem illustrates several options available for 
plotting. Rigid Format 1 - Heat (Static s) was used. The combination of 
parameters used is shewn in the table below. 


Plot Set Axes View Normals Label PTITLE S^ 01 


1 

1 

X Y Z 

Default 


x 

X 

x 1 

2 

2 

X Y Z 

Default 



X 

x 

3 

3 

X Y Z 

Default 

X 

X 

X 

x 

4 

4 

X Y Z 

Default 

• . 


X 

x 

5 

2 

X MYMZ 

Default 




X 

6 

2 

MZ XMY 

0. ,0. ,0. 


X 

X 

x 

7 

2 

Y X MZ 

0. ,0. ,0. 


X 

X 

X 

A sample plot 

is reproduced in 

Figure 6. 


• 




Variance Analysis, Transient 

(DEM03) 



The model 

is exactly the same as was 

used 

for DEM01. 

Initial 

conditions 

were chosen to 

be the 

estimates 

of DEM01 , 

with 

steady loads 

. The 

asymptotic 
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values should agree with the DEM01 results. Runs 3 and 4 show the check- 
point and restart phases. No deviation data was produced. 

A special modeling technique was used to enforce transient boundary 
conditions. The generally accepted method (not used here) is with a good 
conductor to "ground" plus a large heat. The illustrated method uses 
extra points plus DMIG. The "load" applied to the extra point is the 
enforced temperature. The "temperature" or solution of the extra point 
is the heat required to maintain the specified temperature. 


- 17 - 



a 


lo^s’&xna 

ALL. THE ChBQr ElEh^nTS WITH NOHMM.B 

^ 00 



Sample Plot of HBDY Elements 
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